Disorder induced metallicity in amorphous graphene 
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We predict a transition to metallicity when a sufficient amount of disorder is induced in graphene. 
Calculations were performed by means of a first principles stochastic quench method. The resulting 
amorphous graphene can be seen as nanopatches of graphene that are connected by a network of 
disordered small and large carbon rings. The buckling is minimal and we believe that it is a result 
of averaging of counteracting random in-plane stress forces. The linear response conductance is 
obtained by a model theory as function of lattice distortions. Such metallic behaviour is a much 
desired property for functionalisation of graphene to realize a transparent conductor, e.g. suitable 
for touch-screen devices. 



Since the discovery of graphene, and all its unique 
physical properties flPHIl , attention is now turned to func- 
tionalization of this material to suit specific applications. 
One example is the suggestion of graphane, where H 
atoms are adsorbed on graphene. Graphane, which is an 
insulator, was predicted from first principles theory Q, 
and was subsequently realized experimentally Q. Fur- 
ther theoretical works [sl, loj addressed e.g. the value of 
the band gap, which was found to be of order 5.7 eV. 
In a similar fashion, fluorine can also be found to ad- 
sorb, and in theoretical works a band gap of 7.4 eV is 
found fiol [ill , which is larger than the measured value 
of 3.4 eV The reason that a large band gap opens 
up when hydrogen or fluorine is adsorbed on graphene is 
that sp^ bonded C atoms become sp^ bonded. Hence it 
seems possible to decrease the conductivity by chemical 
functionalization, and turn the semi-metallic graphene to 
a semi-conductor or even an insulator. Smaller values of 
the gap can be reached by a functionalization by organic 
molecules [isj . 

The ways to increase the conductivity by chemical 
means has also been discussed, although here signifi- 
cantly less success can be identified. Adsorption by sin- 
gle impurities [14] have been explored, as well as the re- 
placement of C atoms for B or N atoms [l5| , or even the 
creation of structural defects in the C matrix Til- isj]. 
The electronic structure of non-crystalline graphene, e.g. 
around grain boundaries have also been under consider- 
ation, both from tight-binding analysis 1 1911 and within 
the framework of first principles theory [20[ in combi- 
nation with calculations of transport properties. In the 
works mentioned above, the conducting properties are 
still those of a doped semiconductor, with relatively few 
charge carriers. A fully metallic behavior was however 



suggested in an amorphous structure of graphene [2l|, 
where Stone- Wales defects were introduced into graphene 
and geometry optimization was done according to a 
Keating-like potential. After geometry optimization, the 
electronic structure was calculated from a tight-binding 



Hamiltonian, and a non-zero value of the density of states 
(DOS) at the Fermi level (Ep) was obtained, suggest- 
ing that metallic conductivity is possible. The possibil- 
ity of a transparent, mono-atomic thin material, has a 
great potential in applications involving touch screens 
and electronics [22]. Very recently, the two-dimensional 
amorphous carbon has been derived experimentally by 
electron beam irradiation [23]. In this paper we study 
its electronic structure with a first-principles based ap- 
proach, which accurately described both the chemical 
bonding and structural properties, as well as the elec- 
tron energy spectrum. 

The amorphous graphene structures were obtained by 
means of a stochastic quenching method 24 , [2^ . The 
method is well adapted to describe amorphous structures 
that have been obtained by ultra-fast cooling, such as e.g. 
sputtered amorphous films where the structure has been 
realized directly from the gas phase [2^. It can also be 
used to describe the most typical amorphous structures 
of a material in order to find reliable parameters for a 
Hamiltonian description of the liquid state [27']. Nor- 
mally, this technique is used for calculations of amor- 
phous structures in bulk where the atoms are first placed 
randomly in the calculation cell and then the positions 
are relaxed by means of a conjugate gradient method un- 
til the force on every atom is zero. In the present case we 
confined the atoms to the plane in order to find a reliable 
description of the two-dimensional amorphous structure. 
First, 200 atoms were placed at random positions in the 
plane with an average density equal to that of graphene. 
We then relaxed the positions of all atoms, while enforc- 
ing that all atoms should stay in-plane. Then the area of 
the plane was relaxed until the pressure upon the cell was 
zero. At this point we obtained the structure called "ba- 
sic planar" . The area relaxation resulted in new forces 
on the atoms so we had to repeat the position and area 
relaxations until both the pressure and forces were min- 
imized. At this point we obtained the structure called 
"final planar". The relaxation of positions was then al- 
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lowed to take place in 3 dimensions and the process of 
area and position relaxation cycling was repeated until 
the in-plane stress and all forces on the atoms were min- 
imized again. At this point we obtained the structure 
called "buckled". 

The first principles, self-consistent electronic structure 
calculations were performed by means of density func- 
tional theory Q using the projected augmented wave 
method [i^ as implemented in the Vienna ab initio sim- 
ulation package (VASP) code [s^- The geometry opti- 
mizations have been performed without any symmetry 
constraint and using the conjugated gradient algorithm. 
A plane wave energy cutoff of 300 eV was used for the 
structural relaxations whereas 500 eV was used for the 
DOS and total energy calculations. The Perdew-Wang 
3l| parametrization of the generalized gradient approxi- 
mation of the exchange-correlation interaction was used. 
For sampling the irreducible Brillouin zone the F point 
was used for the relaxations whereas a 5x5x1 mesh of 13 
k-points was used for the DOS and total energy calcula- 
tions. The amount of vacuum in the supercell was chosen 
to avoid artificial surface-surface interactions (24 A). In 
all relaxations, convergence was achieved when the force 
on each ion was less than 10^^ eV/A. 

The basic planar configuration that was obtained after 
the first position-area relaxation cycle is shown in Fig- 
ure [TJa) . We can see that this structure consists of 23 
hexagonal and 56 pentagonal, some tetragonal and trian- 
gular rings. There are also some areas with larger rings 
of up to 8 atoms. This structure is still very compressed 
since the area is still close to that of graphene. We can 
see that a large number of carbon atoms are both 4 and 
5-fold coordinated, which indicates a frustrated state of 
high energy. The total energy per atom of this structure 
is 1.50 eV larger than the graphene reference. Further it- 
erations of positions and area relaxations resulted in the 
final planar structure shown in Figure [TJb). The number 
of hexagonal rings has now increased to 38 and the num- 
ber of pentagonal rings decreased to 25. Several larger 
rings have appeared with up to 9 atoms. All carbon 
atoms are now 3-fold coordinated and hence sp^ bonded, 
except in 2 cases where the coordination is 2-fold with 
a 180° bond angle and hence the bonding is sp. There 
is also one atom in the lower edge of the figure that is 
bonded to only one other atom. The total energy per 
atom of this structure is now 0.50 eV larger than the 
graphene reference. 

The position and area relaxation of the final planar 
structure was taken further by letting the atoms relax 
also in the direction normal to the plane. The resulting 
buckled structure is shown in Figure [UJc) and a side view 
can be seen in Figure [Ijd). From the side view we can see 
that the plane has evolved into a weakly oscillating shape. 
In this structure the number of pentagonal and hexago- 
nal rings has not changed compared to the final planar 
structure. The atom with only one bond is also present in 
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FIG. 1: (Color online) (a) Top view of the basic planar struc- 
ture. Hexagonal rings are marked with blue and pentagonal 
rings are marked red. (b) Top view of the final planar geom- 
etry, (c) Top view of the buckled geometry, (d) Side view of 
the buckled geometry. The distance from left to right in (d) 
is 41 A and the buckling amplitude is about 1.7 A. 



this structure and it can be seen in Figure [IJd) where it 
protrudes out below the surface. The total energy of the 
buckled system was 0.43 eV larger per atom than the to- 
tal energy of the graphene reference. This structure must 
be considered as an archetype structure amonst many 
possible realizations of amorphous graphene. There may 
even exist lower energy amorphous states than the one 
we have found here. The amorphous state that we have 
found is surprisingly flat and contains both squares and 
triangles. We speculate that the flatness is due to an av- 
eraging of defect-induced buckling forces. Square geome- 
tries have recently been observed in defect rich graphene 
that was created by electron irradiation [23| . 

The radial and angle distribution functions are shown 



in Figure 2(a) We can see that the first peak in all radial 
distribution functions is located at 1.4 A. However, the 
radial distribution function of the basic planar structure 
has somewhat broader peaks than the other structures. 
The angle distribution functions show a much larger dif- 
ference. Here the basic planar structure shows three pro- 
nounced peaks around 60, 110, and 160 degrees whereas 
the final planar and buckled structures show one pro- 
nounced peak around 120 degrees and only small quan- 
tities of 60 and 90 degree angles. 

The total DOS for graphene and the three amorphous 
geometries are shown in Figure 2(b) We can see that 



the basic planar structure has somewhat broader bands 
but the widths of the final planar and buckled struc- 
ture are somewhat narrower than graphene. The largest 
difference is however around the Fermi level where the 
amorphous structures have high enough density of states 
to be clearly metallic. 
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FIG. 2: (Color online) (a) The radial distribution functions 
(RDF) for all structures together with the angle distribution 
functions (ADF). (b) Density of states of the three amorphous 
geometries compared to the graphene structure. 



The reahzation of a metalhc DOS for a monolayer- 
thick graphene material is a significant result, and the 
transport properties of this system requires further anal- 
ysis. We do this by a model theory, as described below. 
We consider the case of an imperfect 2D carbon lattice 
where all C atoms are interconnected to three other C 
atoms and consider it to be bipartite. In the spirit of a 
nearest neighbor model, we have as a Hamiltonian of the 
system 
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H = V^(k,k')a],fek' +i?.c. 



kk 



(1) 



where ^(k, k') = -{t/N) J2{ij) e^**'''"'+*'''''"^ operators a 
and b annihilate electrons in sublattices A and B, respec- 
tively. At each site i, we express the spatial vectors to the 
neighboring sites rj , can be expressed in terms of and 
additional piece Sij, i.e. Vj = Vi + dij. Assuming that the 
random vector Sij connecting sites i and j are indepen- 
dent of the specific location in the lattice, we can simplify 
the random vector to depend only on j. This allows the 
calculation ^(k,k') = -it/N)J2{^j)e-'^^-^"^■'''+'^'■^' = 
—tS{k — k') 2 3 e^^'^' , so that we can define 0(k) = 

(5(k — k')^(k, k'). The Hamiltonian is thereby reduced to 

H =^0(k)aj,5k + i/.c. (2) 

k 

Next, we assume that the random vector Sj, j — 1,2,3, 



is normal distributed with mean m, 



{mjx,mjy) and 



variance cTj = {ajx,<Jjy)- Also, by assuming that the 
random variables 5jx and 5jy are independent, we can 
replace the exponential e^^'^' by e* '"^ "'^3-''-/^"'^?«''«/^. 
As mean values of the random vectors, we take a set of 
basic lattice vectors for graphene, e.g. mi = a(-v/3, l)/2. 



m.2 = a(— \/3, l)/2, and ma = a(0, — 1), and we assume 
that Ujx = (Jjy = (Jj . For structures with only small vari- 
ations from the perfect graphene lattice, we can expand 
around the K points ±K = ±27r(-\/3/3, l)/3a, around 
which the dispersion relation for graphene is linear, i.e. 
£;g,(k) = |(/)gr(k±K)| = |±t)Ffce±*('^/3-¥')| ^ ypk, where 
vf = iat/i and t&rup = ky/kx- Expanding the potential 
0(k) around ±K, retaining contributions of quadratic or- 
der in k at most, we obtain 



0(k ± K) « ± {i/c 1 =F e 



^±i(7r/3-ip) 



(3) 



where kr 



87r/3A/3( 



where we have introduced a — (27r-\/3cr^/9a) cos(95 — 7r/3) 
and V = iii? exp{— 6(27r(T/9a)^}, whereas cr^ = X^j'^'j- 
The scattering potential converges to the usual graphene 
potential as the variance ct^ — > 0, as required. The 
above expression for the scattering potential holds for 
small distortions from the ideal graphene lattice. Quan- 
titatively we require that a <^ Uu — \J 3\/3a/47rfcc, 

is a high energy cut-off, giving 

cr„ = 3a^y3/27r3/2 ~ 0.4a. 

Numerically, we find that the dispersion given in Eq. 
dS]) is applicable for lattice distortions up to about cr/a < 
0.5, see Fig. [3] (a), where we plot the computed bands 
when applying the stochastic arguments to the model 
given in Eq. ©. Around K, the dispersion is basi- 
cally linear for rather large lattice distortions, and the 
DOS resulting from the band structure. Fig. [3](b), shows 
the typical linear characteristics around the Fermi level 
{Ep = 0) for small lattice distortions {a /a ^ 0.5). Here, 
the DOS (cx ImG"") is calculated in terms of the Green 
function G""(k,a;) ^ N-i[^t(k)(^(k)/(a; - £(k) -|- ir])] is 
defined in terms of the eigenbasis {ip {]<.),£{]<.)} given by 
the Hamiltonian in Eq. ©. In Fig. [21(c), we plot the 
DOS resulting from the approximate dispersion in Eq. 
([3]), showing a monotonic increase with the lattice dis- 
tortion, as expected. 

For further investigations of the potential metallicity 
of the amporphous graphene, we also calculate the linear 
response conductance, Gq, using the Kubo formalism. 
For an electric field applied in the x-direction, we use 



1 /"^ 

Go = lim - / ([jt(q,r), j,(q,t)])e'-^cir, 



(4) 



where limq^o jx(q, t) = «e *(k)aj^6k + H.c, and 



2,3 



xe"" . Neglecting vertex cor- 
rections, we calculate the d.c. conductance at zero tem- 
perature according to 



Go — 



J |$(k) ImG''(k,cj = 0)|2dk/(27r)2. 



(5) 



We obtain a monotonically increasing conductance as the 
distortion of the lattice increases, see Fig. |3] (d), which 
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FIG. 3: (a), (b) Band diagram and DOS calculated within 
the stochastic model based on Eq. for varying degree of 
lattice distortion, (c) DOS{Ef) calculated using the approx- 
imate dispersion relation in Eq. and (d) linear response 
conductance, at low temperatures, based on the calculations 
in (a) and (b). The residual DOS and conductance at zero 
lattice distortion in panels (c) and (d) is due to an artificial 
finite broadening {rj = 0.05) introduced in the Green function. 



supports the conclusion that our two-dimensional amor- 
phous carbon lattice is metallic. The same conclusion is 
drawn from the approximate model based on the disper- 
sion relation in Eq. ([3]). 

In summary, we have performed first principles theory 
as well as considered a model Hamiltonian, to demon- 
strate that graphene can become metallic, once a signifi- 
cant disorder is introduced in the material. In the present 
work, this has been achieved by a realistic amorphization 
in 3 dimensions, using the stochastic quenching method. 
Our results show a transition from the complex electronic 
structure associated with perfectly crystalline graphene 
to a regular conducting behavior. This transition occurs 
rather drastically as a function of the amount of disor- 
der in the materials. Experimentally, disorder and non- 
crystallinity is commonly obtained by ion-irradiation of a 
crystalline material, or by rapid thin-film growth, and it 
is likely that this is an experimental way forward also in 
this case. The energy associated with some of the amor- 
phous structures considered here is only 0.43 eV/atom 
higher than that of crystalline graphene. 

If realized experimentally, a much desired metallic and 
transparent thin film has been found, with the possibil- 
ity to have impact in applications, e.g. in touch screens. 
Given that other allotropes of C can become supercon- 
ducting, at least when doped, this may also be explored 
for amorphous graphene. 

Calculations were performed on Ainil, the supercom- 
puter at the physics institute at UACH funded by 
Chilean FONDECYT projects 1110602 and 11080259. 
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